Density-functional study of adsorption of isocyanides on the gold (111) surface. 
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Density functional theory within the generalized-gradient approximation is used to study the 
adsorption of the isocyanides CNH and CNCIfa on the gold (111) surface at several coverages. It 
is found that these molecules are highly selective in their adsorption site preference. Adsorption is 
possible only at the top site, although binding is rather weak and dipole-dipole repulsion prevents 
binding at coverages of 1/3 ML and higher. At all other high symmetry sites considered (hep, fee, 
bridge), the isonitriles are not bound. To reveal the main mechanisms of bonding, the isonitriles are 
compared to CO and ammonia using Au-X radicals as model systems. The systematic trends are 
understood in terms of a donation from the ligand lone pair and tt* back donation. Finally, CNH 
is found to be more strongly bound to an undercoordinated gold atom (adatom) on the surface. 



I. INTRODUCTION 



Gold is widely chosen as an electrode material for stud- 
ies of molecular conductance for a number of practical 
reasons, including ease of deposition and patterning and 
presentation of a surface that is relatively easy to clean 
and maintain. Another important factor has been the 
extensive base of experience in forming self assembled 
monolayers (SAMs) of organic molecules, linked to the 
gold surface via a sulfur atom. The initial discovery of 
the thiol-gold route to self assembly^ led to rapid growth 
in the study of the formation and properties of organic 
layers^''^. The S-Au bond is relatively strong (~1.7 eV), 
but the barrier to lateral motion between adsorption sites 
is rather low (~0.1 eV), facilitating formation of ordered 
monolayers^. Sulfur has been alsmost exclusively used 
as the "alligator clip"** in the experiments probing con- 
ductance at the molecular scale^'^. However, the sub- 
stantial variation in the measured conductance for ba- 
sic molecules such as alkanes contacted using the Au-S 
link, between different groups and different measurement 
techniques, remains a puzzle. One possibility is sample 
to sample variation in the local structure of the Au-S 
link^. This is consistent with the low barriers for mo- 
tion that facilitate assembly and the known mobility of 
Au atoms during formation of the thiol linked la.yev^. 
Also, the atomic scale structure of the interface of as- 
sembled layers with Au-S linkages has proved to be a 
very delicate problem for ab-initio methods^, although 
a recent study appears to reconcile extensive data for 
decanethiok-''^ . Overall, despite wide usage, the Au-S 
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link displays a number of critical issues for application in 
molecular conductance studies. Recent research has been 
directed to find alternative link chemistries, e.g. Ru or 
Mo carbeneiiii?. 

An alternative link chemistry for a gold electrode that 
shows more selectivity in the local bonding geometry 
is an attractive possibility in the context of molecu- 
lar scale conductance. Although Au surfaces are rela- 
tively unreactive, other candidates have been explored. A 
very recent study demonstrated the utility of amine ter- 
minated molecules for molecular conductance studiesi^. 
The AU-NH2-R links proved to be flexible and repro- 
ducible. Another example, with a similar bonding char- 
acter, is isocyanide (CN-R, e.g. CNH), a less-stable 
isomer of cyanide (NC-R, e.g. NCH). Well-established 
synthetic methods exist for the CN- linker with vari- 
ous molecules (R). Bonding to Pt or Au electrodes has 
been observedi^. In contrast to the Au-S link, adsorp- 
tion of isocyanide (Au-CN) linked molecules on gold sur- 
faces is much less studied. Most experimental works 
on this subject consider adsorption of the molecules on 
gold nanoparticlesi^4iS4i2ii&. Henderson et alM- stud- 
ied the formation of self-assembled monolayers of diiso- 
cyanides on the gold (111) surface. It was concluded 
from reflection absorption infrared (RAIR) spectra that 
the molecules with rigid chains between two CN- groups 
(such as biphenyldiisocyanide) attach to the gold sur- 
face through only one isocyanide group while the other 
group remains free. EUipsometry data indicate that in 
these cases, the molecular axis is perpendicular to the 
surface. Molecules with flexible alkyl chains can attach to 
the gold surface through both isocyanide groups, which 
makes them unlikely candidates for molecular conduc- 
tance studies. The preferred adsorption site and geome- 
try remain as unresolved questions in the literature. An- 
gelici et al. observed the peak in the infrared spectra 
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associated with N-C stretch mode to be shifted 50-70 
cm~^ higher upon adsorption on gold powderii^ii^. They 
concluded that isocyanides bind to the surface at the top 
site. Henderson et al}^ observed a similar blue shift for 
adsorption on the gold (111) surface, but concluded that 
both atop and three-fold adsorption configurations are 
possible. 



The earliest of the few theoretical works on the ad- 
sorption of isocyanide on gold considered phenyl iso- 
cyanide attached to a single gold atom, a gold dimer 
and a trimep^o. The calculations showed that for each 
gold cluster, the binding was to a single Au atom. Bond- 
ing motifs analogous to the bridge and hollow sites on 
the (111) surface were not observed. The results sug- 
gest site selectivity for the isocyanide link, although the 
details of the adsorption energy and the geometry on 
the Au(lll) surface remain to be determined. Recently 
several groupsSliSS successfully fabricated self-assembled 
monolayers of various diisocyanide molecules sandwiched 
between two gold electrodes, and studied electrical trans- 
port through them. The impact of the isocyanide link on 
the conduction has received some attention from theorists 
as well. Conductance and voltage-current characteristics 
of molecules connected to metallic electrodes were cal- 
culated using combined density functional theory (DFT) 
and non-equilibrium Green's function methodsS^. How- 
ever, the adsorption geometry was not explicitly deter- 
mined by energy minimization. 



This paper presents a study of the adsorption of iso- 
cyanide linked molecules on the gold (111) surface using 
a density functional theory (DFT) approach. We focus 
on the binding of the essential link element to the metal 
surface, considering the simplest CNH molecule. 

Calculations performed for CNCH3 confirm the results. 
The basic electronic structure of the CNH molecule sug- 
gests that the filled lone pair on the C atom will be avail- 
able for donor-acceptor type bonding to Au complexes, 
in analogy with ammonia (NH3) and CO. We chose to 
draw a comparison between CNH, CO and NH3 in order 
to better understand the factors that affect adsorption. 
The case of CO has been extensively studied for many 
transition metal surfaces24. There are also recent calcu- 
lations for the case of NH3 on Au(lll)S^. The standard 
Blyholder picture2& describes the bonding in terms of the 
trade-off between a donation from the lone pair to the 
metal and back donation from the metal to the empty 
TT* molecular orbital. Two main themes emerge: (1) the 
interplay of the electronic couplings in the Blyholder pic- 
ture with the corresponding charge transfers; (2) the im- 
pact of the dipole moment, both that of the molecule 
prior to adsorption and the induced dipole upon adsorp- 
tion. In the next section, the methodology and the basic 
results of the DFT calculations are described. The phys- 
ical interpretation of the results is discussed in Sect. HI. 



II. METHODOLOGY AND RESULTS 

The adsorption of the target molecules on the Au sur- 
face is studied using a computational approach based on 
DFT. The main calculations are done using the general- 
ized gradient approximation (GGA) of Perdew, Burke 
and Ernzerhof (PBE)22.. Different forms of GGA are 
widely used for metal surface and adsorption problems 
in order to address the large errors in adsorption en- 
ergy found with the local density approximation (LDA), 
although the accuracy still varies with the system un- 
der studySiiSSiSSi^fli^. For comparison, selected calcula- 
tions are also done with the LDA224^. To have an ini- 
tial picture of the bonding, the target molecules bonded 
to a single Au atom are considered using the NRLMOL 
DFT codc^'^. For the light atoms (C, N, and H), an all 
electron basis is used. A slight modification'^^^ of norm- 
conserving TrouUier-Martins pseudopotentials is used for 
gold atoms. NRLMOL uses Gaussian basis sets opti- 
mized for density-functional calculations^. The expo- 
nents and contraction coefficients are determined by op- 
timizing the total energy of the free atom in the ground 
state. For the valence states, the number of independent 
contractions used (s, p, d) for each atom is as follows: H 
(4, 3, 1), C (5, 4, 3), N (5, 4, 3) and Au (4, 2, 4). It has 
been shown that these basis sets are well-converged and 
that they have negligible basis set superposition error. 
Further details of the basis sets can be found in Ref.— . 
Radicals were treated with spin-polarized, unrestricted 
calculations. 

The adsorption on the surface was modeled with the 
standard scheme of a periodically repeated slab geom- 
etry. We found that the slab of four layers of gold 
atoms was enough to reproduce key adsorption charac- 
teristics such as the binding energy. The molecules are 
adsorbed on one side of the slab. For most of the cal- 
culations, the molecular structure is relaxed, with the 
molecular axis kept perpendicular to the surface and the 
Au atoms in the slab frozen at their bulk positions. Re- 
laxation is carried out until the maximum force is less 
than 0.15 eV/A. For selected cases, the top layer of Au 
atoms was allowed to relax. We also explored bonding 
of CNH to an Au adatom in the hep hollow site, al- 
lowing the adatom and the surface Au atoms to relax. 
Most of these calculations were done using the WIEN2k 
DFT code (full potential and linearized augmented plane 
waves'^'*). For comparison, some calculations were also 
done with the ABINIT package?^''*^, which employs a 
plane wave basis set and pseudopotentials. LDA cal- 
culations were done with Hartwigsen-Goedecker-Hutter 
(HGH) pseudopotentials^, while GGA calculations were 
done with TrouUier-Martins type pseudopotentials^Si^. 
It was checked that the energy is converged with respect 
to the number of k-points as well as the size of the basis 
set. To facilitate analysis, the partial densities of state 
(PDOS) were plotted in selected cases. In the WIEN2K 
code, a modified tetrahedron method'*^ is used for density 
of state calculations; the PDOSs are calculated by pro- 
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jecting the wave functions onto the spherical harmonic 
basis functions centered on the atoms in question inside 
a sphere centered on the atom (atomic sphere radii are 
1.06 A for Au atoms, 0.58 A for C and N, and 0.29 A for 
H). 

We start by considering binding of the CNH and 
CNCH3 molecules to a single gold atom. The binding en- 
ergy of the Au atom, the gold to ligand distances, and the 
dipole moments of the free molecules and the molecules 
bound to the gold atom are reported in Table I. The re- 
sults for CNH and CNCH3 are very similar. For compar- 
ison, results for CO and NH3 are also shown. The first 
step of relaxation allows only configurations with axial 
symmetry. The axial CO and CNH molecules are very 
similar in their electronic structure. The highest occu- 
pied molecular orbital (HOMO) is the lone pair on the C 
and there is a doubly degenerate tt bond between C and 
the N or O. For NH3, the N lone pair is the HOMO, but 
there are no tt bonds. The a part of the interaction in 
all three is due to coupling of the lone pair orbital on the 
ligand to the d^2_^2 and s states on the Au. This results 
in partial transfer of electron density onto the Au atom. 
Because the Au s state is half filled, the frontier, anti- 
bonding a state in the complex is also half filled. The 
initial ligand lone pair orbital energy is highest for NH3, 
comparable to the Au s state, and successively lower for 
CNCH3, CNH and lowest for CO. 

The a donation is partially balanced by the back dona- 
tion from the Au dzx and dyz into the empty tt* doublet 
on the molecule for CO, CNH and CNCH3. This balance 
is slightly different for AuCNH and AUCNCH3 relative to 
AuCO. In the AuCO case, the C-0 bond length increases 
by 0.004 A relative to CO. In the AuCNH (AUCNCH3) 
case, the C-N bond length decreases by 0.007 A relative 
to CNH (CNCH3). This difference suggests less tt* back 
donation in the AuCNH case. The reduced back dona- 
tion is consistent with the tt* doublet being about 1 eV 
higher in CNH and even a bit higher for CNCH3. This 
is also consistent with the trend of increasing Au-C bond 
length from AuCO to AuCNH to CNCH3. The trend for 
net charge transfer to the Au as measured by MuUiken 
type analysis is consistent with larger a donation and 
smaller tt back donation from AuCNCHs to AuCNH to 
AuCO. The binding energy follows the same order, but 
the impact of larger a donation and smaller tt* back do- 
nation are competitive, with the cr donation dominating 
the trend. By comparison, the NH3 molecule shows the 
strongest a donation, with the largest charge transfer to 
Au, but offers no empty tt space for back donation. Its 
binding energy to a gold atom is 0.60 eV, less than for 
CO. Correspondingly the Au-N bond length is substan- 
tially larger (2.277 A). The role of the tt states is further 
highlighted by the results of full relaxation, which leads 
to bent structures for Au-CNH, AU-CNCH3 and Au-CO. 
This bending is due to the pseudo Jahn- Teller effect. The 
lowered symmetry allows mixing between the half occu- 
pied cr antibonding level in the complexes and one of the 
empty tt* states. The resultant energy gain depends on 



the gap between a and tt* states. It is largest for the Au- 
CO complex and smallest for AU-CNCH3; the Au-CO 
complex has the smallest gap between tr and tt* states, 
while AU-CNCH3 has the largest. 

In comparison to CO, a key difference is that CNH and 
NH3 have a much larger dipole moment. Table 1 shows 
the calculated dipole moments which agree well with the 
measured values. In all three cases, there is a substantial 
increase in the dipole of the gold-ligand complex relative 
to the isolated ligand. The dipole moment of the AuCO 
complex is 1.6 D for the linear configuration, compared 
to the small dipole of the CO molecule. The dipole of 
the AuCNH complex increases by 2.4 D over the CNH 
ligand, while the increase is 3.2 D for the NH3 case. The 
increase in dipole can not be simply seen in a point charge 
approximation with associated charge transfer (e.g. from 
a MuUiken type analysis). Instead, the extension of the 
lone pair electrons at the C end of the ligand onto the 
Au by quantum mechanical mixing with the Au s state 
is responsible, partly balanced by the tt* back donation. 
The reduced back donation in the AuCNH case results 
in a larger change in dipole. In the NH3 case there is 
no back donation and the dipole change is even larger. 
On the other hand, the distortion to the bent form in 
the Au-CNH and Au-CO cases results in a much reduced 
dipole, consistent with the enhanced back donation from 
the Au to the C. 

In addition to the hybridization considerations, the 
dipole also infiuences the relative binding of Au to the 
ligands. There is an additional attraction in the Au-CNH 
complex due to polarization of the Au by the dipole. This 
influences the differences in binding energy, in addition 
to the systematic hybridization differences. To estimate 
this contribution, we use the calculated polarizability of 
the Au atom.^- together with a point charge model of the 
molecule derived from fitting to the self consistent elec- 
trostatic potential and the calculated dipole moment. For 
the case of AuCNH, we estimate that polarization con- 
tributes 0.1 eV to the binding energy. For AuNHs, the 
effect is even smaller (0.01 eV). 

So far as we know, there are no direct measurements 
of the Au binding energy in these radicals. The Au- 
NH3 radical was studied using both DFT approaches and 
quantum chemistry approaches for evaluation of the cor- 
relation energy (e.g. coupled cluster with singles, doubles 
and triple corrections, CCSD(T))^^ . Substantial care was 
also taken to account for basis set superposition errors, 
which can be a significant issue for these correlation tech- 
niques in particular"*^. The Au binding energy in this 
case was deduced to be 0.78±0.06 eV, similar to the 
GGA-PBE value we obtain. There are several studies 
of the Au-CO radical with DFT and quantum chemistry 
approachea^Si^SiiSiffi. Our result for the lowest energy ge- 
ometry, the bent form, agrees reasonably well with pre- 
vious values based on DFT with CCA: 0.71 eVi^ and 
0.80 eV^°. However, the correlated electron calculations 
consistently give smaller binding energies, ranging from 
0.14 to 0.56 eV depending on the detail a^^i^^ . Given con- 
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Distance, Ang 

FIG. 1: Total energy (GGA-PBE) per unit cell for HNC 
monolayers with 1/1 and 1/3 coverages on the gold slab vs. 
distance d to the surface. Molecules are at the top or filled hol- 
low site with molecular axis perpendicular to the surface. The 
indicated adsorption energy is relative to a uniform molecular 
layer well separated from the surface. 

cerns about basis set superposition errors, it is premature 
to draw a final conclusion from tlie quantum chemistry 
calculations in this case. 

Turning to the adsorption on the surface, we initially 
consider a full monolayer of CNH adsorbed on one side 
of a slab of four gold layers. The supercell contains one 
molecule and four gold atoms, two of them being surface 
atoms. First, the molecules are placed at the top site 
(above the surface Au atom) and the distance from the 
carbon atom to the surface is varied. The molecular axis 
is kept perpendicular to the surface. The observed de- 
pendence of total energy vs. distance d is shown in Fig. 1. 
The binding energy Eads is the energy per molecule re- 
quired to separate the layer of CNH from the Au (111) 
surface, leaving the free standing CNH layer otherwise 
unchanged. The value of Eads is found to be 0.43 eV and 
the optimal distance is 2.01 A. In order to check that 
four layers of gold are enough to represent the surface, 
we repeated the calculation for three, five, and six gold 
layers in the slab. The results value of Eads varies from 
0.453 to 0.431 eV, indicating that the choice of four layers 
is justified. 

When a molecule is placed at the filled (hep) hollow 
site, the energy vs. distance (see Fig. 1) is qualitatively 
different from that at the top site. Although there is a 
local minimum at around 1.6 A, the energy is minimal 
when the molecule is far away from the surface. This 
indicates that the CNH molecule cannot be stabilized at 
the hollow hep site. The same result is obtained for the 
empty hollow (fee) and for the bridge sites. 



The next candidate for a stable configuration at uni- 
form commensurate coverage is a monolayer with one 
molecule per three surface atoms, forming a x 
^/3R30° equilateral triangular lattice on the surface, 
where R=2.883 A is the distance between nearest neigh- 
bor gold atoms. The distance between molecules in this 
case is 4.944 A. For the CNH molecules positioned at the 
top site, we relax the monolayer keeping the molecular 
axis normal to the surface and gold atoms frozen. The 
adsorption energy is found to be 0.18 eV relative to the 
separated slab-monolayer system. For the hep site and 
\/3 X \/3R30° surface unit cell, similar to the 1x1 case, 
the energy is minimal when molecules are far away from 
the surface, although there is a local minimum in energy 
curve (see Fig. 1). This shows that adsorption is not 
possible at the hep site. Similar results are obtained for 
the CNH monolayer with 1 /4 coverage (2x2 surface unit 
cell) and the CNCH3 monolayer with 1/3 coverage (see 
Table H). 

So far binding of pre-formed monolayers was consid- 
ered. For adsorption from the gas phase, the energy cost 
of forming the free monolayer must be taken into account. 
The CNH molecule has a large dipole moment, as noted 
above. In the case of 1 ML coverage the distance be- 
tween molecules in the monolayer is 2.883 A. The point 
dipole approximation gives a crude estimate of about 1.4 
eV per molecule for the energy cost of monolayer forma- 
tion. This is much larger than the energy gained upon 
adsorption. Therefore HNC molecules cannot adsorb on 
the (111) gold surface at 1 ML coverage. For smaller cov- 
erages (1/3 and 1/4 ML) the distance between molecules 
is large enough for the point dipole approximation to be 
accurate. We calculate an energy cost of 0.26 eV per 
molecule to form a CNH monolayer at 1/3 ML coverage 
and 0.17 eV at 1/4 ML coverage. When compared to 
the adsorption energies in Table II, the formation of 1 /3 
ML is endothermic, while the 1/4 ML coverage is slightly 
exothermic, about 0.03 eV per molecule. 

Just as noted for the AuCNH complex, adsorption 
leads to a change in dipole moment, due to both the 
charge transfer and the polarization of the molecules and 
the gold surface. For sparse coverages (less than 1/3), the 
induced dipole should not depend much on the distance 
between molecules. This extra dipole gives a coverage de- 
pendent contribution which is included in the adsorption 
energy reported in Table II. The scaling of this term is 
inversely proportional to the cube of the intermolecular 
distance. Using this approach and the data in Table II, 
we extrapolate the adsorption energy in the dilute limit 
to be 0.23 eV and the effective dipole moment of the ad- 
sorbed monolayer to be 3.35 D. Thus the induced dipole 
moment in the monolayer is about 0.3 D, a considerably 
smaller change than found for the AuCNH complex. We 
also estimate the contribution to the binding that derives 
from the classical polarization response of the metal sur- 
face to the initial dipole on the CNH. We use the point 
charge model mentioned previously for CNH and assume 
a classical image potential to arrive at about 0.1 eV. In 
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TABLE L Comparison of Au-X complexes, where X is CO, NH3, CNH and CNCH3. All eharacteristics are calculated with 
NRLMOL DFT code using both PBE or LDA exchange-correlation potentials. Experimental values of the dipole moments are 
given in parentheses when available. The X dipole column corresponds to the dipole moment of the isolated molecule X in the 
case of linear configurations. In the case of a bent configuration, it corresponds to the dipole moment of the ligand X with the 
geomtery taken from the fully relaxed bent Au-X complex. 

Molecule X Au-X geometry Functional Binding energy, eV Au-X bond length, A X dipole, D Au-X dipole, D 



CO 


linear 


PBE 


0.68 


1.997 


0.20(0.12)" 


1.65 


CO 


linear 


LDA 


1.27 


1.948 


0.23 


1.57 


CNH 


linear 


PBE 


0.93 


2.002 


3.09(3.05)'' 


5.46 


CNH 


linear 


LDA 


1.54 


1.957 


3.20 


5.50 


NH3 


linear 


PBE 


0.60 


2.277 


1.53 


4.72 


NHa 


linear 


LDA 


1.07 


2.188 


1.54 


4.89 


CNCHa 


linear 


PBE 


0.99 


2.011 


3.95(3.83)^^ 


6.97 


CO 


bent 


PBE 


0.87 


2.005 


0.14 


0.35 


CNH 


bent 


PBE 


1.04 


1.989 


2.69 


2.89 


CNCHa 


bent 


PBE 


1.05 


2.019 


3.94 


5.74 



"J. S. Muenter, J. Mol. Spectrosc. 55, 490 (1975). 

'G.L. Blackmail, R.D. Brown, P.D. Godrey, and H.I. Gunn, Na- 
ture 261, 395 (1976). 

'S.N. Ghosh, R. Trambarulo, and W. Gordy, J. Chem. Phys. 21, 
308 (1953). 



TABLE II: Stabilization energies from WIEN2K calculations 
for various molecules obtained with GGA-PBE exchange- 
correlation potential. 

Molecule Coverage Site Stabilization energy, eV 



CNH 


1/1 


top 


0.43 


CNH 


1/1 


hep 


no adsorption 


CNH 


1/1 


fee 


no adsorption 


CNH 


1/1 


bridge 


no adsorption 


CNH 


1/3 


top 


0.18 


CNH 


1/3 


hep 


no adsorption 


CNH 


1/4 


top 


0.20 


CNH 


1/4 


hep 


no adsorption 


CNCHa 


1/3 


top 


0.22 


CNCHa 


1/3 


hep 


no adsorption 


CO 


1/3 


top 


0.20 



comparison to the relatively small extrapolated net ad- 
sorption energy, polarization plays a noticeable role. 

Given the extra repulsion from the induced dipole, one 
might have expected the CNH to adopt a bent configu- 
ration on the surface. In the radical, this dramatically 
reduces the dipole, as seen in Table I. This has been in- 
vestigated by considering small distortions of the CNH 
away from the vertical a-top structure at 1 /3 ML cover- 
age (molecular axis was tilted by 10°, while the geometry 
of the molecule was otherwise unchanged). Relaxation 
of this initial configuration results in an essentially ver- 
tical structure. Alternatively, we have started from a 
bent CNH configuration at 1/4 ML coverage (with Au-C 
bond vertical, CNH molecule positioned as in the bent 



Au-CNH complex) and relaxed the structure. We find 
a smooth path back towards the vertical geometry. Evi- 
dently the back donation driving force, characterized fur- 
ther in the next section, is weaker in the adsorbed case 
than in the radical. 

In the neutral radical, the net effect of the donor- 
acceptor bonding is relatively weak due to the partially 
occupied Au s-orbital. For an Au cation, or a complex 
which partially withdraws the Au s-electron, the binding 
can be substantially stronger. Experiments and calcula- 
tions have been analyzed for a series of ligands bonding 
to the gold cation^. The estimated bond energies for 
Au+CO, Au+NHg, and Au+CNCHg are 2.1, 3.1 and 3.1 
eV respectively. These are all substantially larger than 
for the corresponding neutral radicals. Furthermore, it 
is known that isocyanides as well as amines and phos- 
phines passivate gold nanoparticlesi^iSiS. This raises 
the interesting question of the character of CNH bond- 
ing to undercoordinated gold atoms on the surface. To 
probe this question, we calculated the binding of CNH to 
an Au adatom using ABINIT. The adatom was relaxed, 
together with the first full surface layer, in the hep hollow 
site of a 4 ML slab with 1 /4 coverage. Then the CNH was 
added to the system and fully relaxed. The final geome- 
try is close to vertical with a reduced Au-C bond length 
(1.958 1) and shghtly longer C-N bond length (by <0.01 
A) relative to the original a-top geometry. The bind- 
ing energy is increased by 0.9 eV, relative to the a-top 
adsorption site on the flat surface. This binding energy 
refers to an array of CNH molecules in the same unit cell 
relative to the gold surface including the adatoms. The 
formation energy of the adatom is not explicitly included 
since we consider the adatom to be a proxy for underco- 
ordinated gold atoms on the surface of a nanoparticle or 
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other rough gold surfaces. 

Finally, it is well known that the DFT method using 
the LDA overestimates binding energies. However it is 
usually expected that general trends are reproduced cor- 
rectly. Some calculations for the radicals were repeated 
with the LDA exchange-correlation functional (Table 1). 
As expected, LDA systematically overestimates binding 
energies by about 0.6 eV and underestimates the Au- 
ligand bond lengths. Also, some of the slab calculations 
of CNH adsorption on Au(lll) were repeated using LDA. 
Adsorption energies calculated with ABINIT for a 1/4 
ML coverage are 0.95 eV, 1.09 eV and 1.12 eV for the 
top, hep and fee sites respectively. Similar results are 
obtained with WIEN2k code for 1/3 ML coverage: 0.80 
and 1.03 eV for the top and hep sites. There is a qualita- 
tive difference between the GGA-PBE and LDA results, 
the GGA being more realistic where experimental tests 
are available. With GGA the top site is preferred and 
there is no adsorption at the hep or fee sites. With LDA, 
on the contrary, adsorption at the hep site is possible and 
it is energetically favored over adsorption at the top site. 



III. DISCUSSION 

We have found that isocyanides selectively bind to the 
atop site on the flat gold (111) surface, with a small ad- 
sorption energy. The C-N bond in the isocyanide is ver- 
tical. This binding geometry is consistent with experi- 
ments described in the introduction: top site adsorption 
is observe di'i^o and ellipsometry data indicate vertical 
orientation of adsorbed isocyanides^^. At the present 
there are no measurements of binding energies or satura- 
tion coverage of isocyanide molecules on the gold surface 
for direct comparison to our calculated adsorption en- 
ergy. The CO molecule, on the other hand, is very well 
studied. As noted, its electronic structure is similar to 
that of the CNH molecule. We performed calculations 
identical to those described above for CO on the gold 
(111) surface and obtained 0.2 eV for the adsorption en- 
ergy, in good agreement with previous calculations, but 
less than the experimental value 0.4 e\—. 

The bonding of both CO and CNH to transition met- 
als is often discussed in terms of the Blyholder picture2&. 
Binding involves donation from the filled lone pair on the 
C into the partially occupied Au s shell together with 2it* 
back donation into the 27r* molecular orbital due to in- 
teraction with gold dxz and dyz- We probe this picture 
for the bonding of CNH on the flat surface by analyzing 
the angular character of the local density of states for the 
adsorbed molecule and the gold slab. The projected den- 
sity of states (PDOS) for the monolayer adsorbed onto 
a Au(lll) slab at 1/3 coverage in the atop configura- 
tion is plotted in Fig. 2 along with the PDOS for a free 
CNH monolayer and a clean Au(lll) slab. Also shown 
are isosurface plots of selected wavefunctions (P-point) 
corresponding to certain energy regions in the PDOS of 
the Au(lll)-CNH system. 



In the PDOS of the free CNH monolayer several very 
sharp peaks are observed which can be associated with 
the molecular orbitals 4(T, Itt, 5ct and 2ti* (the last be- 
ing unoccupied). After adsorption the Aa level docs not 
change much and there is no disccrnable interaction with 
gold. This is not surprising considering that this level 
lies well below the bottom of the gold valence band. In 
contrast, Itt and Str levels undergo substantial reorgani- 
zation. The 5(7 state strongly interacts with hybridized s 
and dz'i orbitals of gold. A localized state is formed below 
the bottom of the gold valence band. In addition, there 
is a state of anti-bonding character that can be identified 
above the occupied gold d-band, but below the Fermi 
energy (around -1 eV in Fig. 2). These states are very 
similar to the hybrids that forms in the AuCNH radical. 
In the radical, the anti-bonding state is half occupied. As 
a result of the strong a interaction, the Itt level now lies 
above 5(T derived level. The Itt level interacts with gold 
dxz and dyz orbitals. This results in two peaks visible 
in the PDOS at around —7 eV and —5 eV, of bonding 
and anti-bonding character respectively, as can be seen 
from the wavefunction plots. The bonding combination 
appears to form a localized state just below the bottom 
of the Au d-bands, while the anti-bonding combination 
is resonant with the d-bands. Finally, the 2tt* peak from 
the CNH derived state density extends somewhat below 
Fermi energy upon adsorption. The corresponding wave- 
function isosurface plot shows that these occupied states 
are of bonding character between gold p-like states avail- 
able near the Fermi energy and the CNH 2tt* orbital. 
Hence back donation into the 2tt* orbital contributes to 
the bonding. Above the Fermi energy, near the center of 
the resonance, the coupling is weak, but clearly of anti- 
bonding character to gold d^z and dyz orbitals. This 
qualitative picture from the PDOS is very similar to pre- 
vious results for CO adsorbed on the a-top site of gold 
(lll)Si. 

It is hard to quantify this analysis of bonding character 
since the PDOS is calculated using only electron densities 
inside spheres surrounding the atoms. However, another 
metric related to bonding trends is the change in C-N 
bond length upon adsorption on gold. The calculation 
for the Au-CNH complex shows a C-N bond length of 
1.173 A in the straight configuration. Upon full relax- 
ation to the bent configuration, the C-N bond is 1.198 
A, compared to 1.176 A in the free CNH. This net weak- 
ening of the C-N bond upon bending traces to the pseudo 
Jahn- Teller effect driven by increased occupation of the 
antibonding 27r* orbital. On the other hand, upon ad- 
sorption on the Au(lll) surface, the calculated C-N bond 
length is slightly decreased to 1.16 A. This contrasts with 
results for CO adsorption, where the C-0 bond length 
slightly increases upon adsorptiorJi. The decrease in C- 
N bond length suggests slightly less tt back donation for 
CNH adsorption on the surface relative to the radical. 
However, careful study of metal-carbonyl complexes re- 
veals a more complex picture^. The final C-N bond 
length (and the corresponding stretch frequency) reflects 
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the balance between a bond polarization due to the metal 
and TT* back donation effects. This balance is affected by 
electrostatics, particularly in the limit of the positively 
charged carbonyl complex. The positive charge near the 
C influences the heteropolar C-0 a bond, driving it to 
be more symmetrical and stronger. This also drives an 
increase in the C-0 stretch frequency in the carbonyls. 
One would expect similar effects on the C-N bond for the 
CNH case. This is consistent with the observed increase 
of the C-N stretch frequency in isocyanides adsorbed on 
goldi^. 

The initial dipole on the molecule and the induced 
dipole upon adsorption both play a role. The CNH 
molecule has a substantial dipole moment, in contrast 
to the CO molecule. Therefore, in the case of CNH, 
the binding energy contains some contribution from the 
polarization energy due to the interaction between the 
molecular dipole and the gold surface. Our estimates 
indicated a modest effect (0.1 eV), which is however a 
non-trivial fraction of the final adsorption energy for an 
isolated CNH on the flat surface (0.23 eV). Furthermore, 
the net dipole inhibits formation of a dense film on the 
surface. The inferred induced dipole on the metal surface 
is substantially smaller than that found in the isolated 
Au-CNH radical. This is indicative of less net charge 
transfer in comparison to the radical. This is not sur- 
prising from two points of view. First, the work function 
of Au(lll) (5.31 eV— ) is smaller than the chemical po- 
tential of the Au atom (conventionally the average of the 
ionization potential and the electron affinity, 5.76 eV— ). 
This inhibits the sigma donation from the lone pair. Sec- 
ond, the surface Au s-state is substantially involved in 
band formation, having nine nearest neighbors. This in- 
terferes with hybridization to the lone pair, reducing the 
net energy gain for the sigma donation process. 

The role of charge transfer has been controversial in 
the analysis of these weakly adsorbed systems and the 
impact of the hybridization implicit in the Blyholder pic- 
ture has been debatedSiSS*Si^ . Bilic et al. argue that 
the main mechanisms of binding of NH3 to the gold sur- 
face are polarization effects and dispersive interactions, 
not the covalent bonding2S,. In their work on pyridine^, 
they find more evidence for charge transfer, but still ar- 
gue that covalent effects are minimal. A more recent 
study of pyridine binding to Au^^ analyzed the impact 
of Au coordination and suggested a more prominent role 
for hybridization effects. The utility of the Blyholder pic- 
ture becomes more apparent when examining trends in 
the binding, e.g. for bonding to different metalsSi. Al- 
ternatively, it gives a way to rationalize the trends for 
the binding of different molecules to the same metal, as 
described for the Au radicals in Section H. Finally, the 
sigma donation becomes more prominent for the case of 
binding to an Au adatom on the surface. We found a 
large increase in binding energy of CNH compared to the 
flat surface (0.9 eV). A similar increase was recently re- 
ported for NH3 (0.4 eV)i2i. The donor-acceptor binding 
to the adatom is also enhanced due to the slight posi- 



tive charge on the Au adatom, similar to the increased 
binding energy for the Au(I) cation to several ligands^. 

The small calculated value of adsorption energy of iso- 
cyanides on the flat gold surface seems to contradict 
the experimental observation of SAM formationiS. Two 
caveats are due here. First, layer formation takes place 
in a solution. Solvation effects can screen the dipole- 
dipole interactions and facilitate layer growth from nu- 
clei. Also, a layer of larger molecules (e.g. alkanes) may 
be more stable than the binding energy of the CNH link 
moiety would suggest due to attractive inter-molecular 
interactions between the extended molecules in the layer, 
e.g. due to dispersive interactions between neighboring 
alkanes. Second, the presence of undercoordinated Au 
binding sites at step edges on the surface may facilitate 
nucleation of the isonitrile layers. We have found that the 
binding energy to the Au adatom is substantially larger 
than the binding energy to the top site on the flat sur- 
face. This may also account for the binding to powdered 
gold samples and nanoparticles where undercoordinated 
Au sites would be common. 

Assessing the accuracy of DFT calculations for these 
weakly bound systems is difficult. The information for 
the Au radicals was already assessed in Sec. H. In the 
case of Au-NHs, the GGA Au to NH3 binding energy is 
smaller than accurate quantum chemistry calculations by 
about 0.2 eV. Bilic et al. studied the binding of NH3 to 
the gold surface using the GGA with the PW91 exchange- 
correlation functional, which in this particular case gave 
an adsorption energy in quantitative agreement with ex- 
periment (within 0.1 eV)^'^. For the Au-CO radical, the 
binding energy is less well established, as discussed in Sec. 
II. Gajdo et al. studied adsorption of CO on (111) surface 
of various transition and noble metalsS4. In particular, 
they report results for CO adsorption on Au and Ag sur- 
faces obtained with PW91 and a revised form of the PBE 
exchange-correlation potentials which has in other cases 
given more accurate adsorption energieaSS. Adsorption 
energies were very sensitive to the choice of the func- 
tional. Both functionals underestimate adsorption ener- 
gies. The revised PBE leads to endothermic adsorption 
while PW91 underestimates the binding energy by about 
0.2 eV. The essential mechanisms of binding are very sim- 
ilar for CO and CNH molecules, as it can be seen from the 
comparison of PDOS pictures. The experience with CO 
on the flat Au surface suggests that the present DFT cal- 
culations may underestimate the binding energy for the 
CNH case as well. Also, in the CNH case, we have found 
that the LDA overestimates the radical binding energies 
and predicts the wrong binding site on the flat surface. 
The role of dispersion forces in the binding has been de- 
bated. The evidence just summarized points to the GGA 
calculations underestimating the binding energy of CO 
and CNH on the Au surface. This is consistent with 
experience in the case of intermolecular interactions for 
closed shell molecules which are dominated by dispersion 
forces. The GGA approximations are known to underes- 
timate binding in these cases while LDA over estimates 
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